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Abstract 

We study the physical origins of the 0(a s ) and 0(a 2 ) corrections to the c — s current 
in the decay b — > ccs in the threshold region 5 = (Mb — 2m c )/2Mb 1. We obtain the 
corrections which are produced by the Coulomb force between the anti-charm and strange 
quarks. The Coulomb corrections Cpir 2 at 0(a s ) and —Cpir 2 ln 5 at 0(a 2 s ) account for 
300% and 120% of the corresponding terms in the Abelian-type perturbative corrections 
respectively. The differences between the Coulomb and perturbative corrections imply that 
other corrections which have other physical origins exist. We also check that the Wilson 
coefficient for the anomalous dimension of the c — s current reproduces the leading and the 
next-to-leading logarithmic terms in the perturbative corrections. 
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1 Introduction 



Disagreement between the experimental measurements and the theoretical predictions for 
the branching ratio of the semi-leptonic decay of the B meson, BK(B —>■ Xli/i), has been 
reported pp. Because the non-perturbative corrections related to the confinement of the 
participating quarks inside the hadrons are at a level of a few percent [TJ, the branching 
ratio of the semi-leptonic decay is mainly determined by the decay rates of the three decay 
modes of the free b quark decay, the semi-leptonic mode b — *> clui, the non-leptonic mode with 
light quarks b — > cud, and the non-leptonic mode with an extra charm quark b — > ccs. Here, 
the non-leptonic modes contribute to the branching ratio of the semi-leptonic decay through 
the total decay width of the b quark. About 30% enhancement of the 0(a s ) corrections 
to the c — s current in the non-leptonic mode b — > ccs has been obtained 0, El IU [£>]■ 
When the large enhancement at the 0(a s ) level and the evolution of the renormalization 
group from the scale of the W boson mass to that of the b quark mass are taken into 
account, the significant disagreement with experiments vanishes [31 HI El [7j. It has also been 
reported that the 0(a 2 ) corrections to the non-leptonic mode b —>■ cud remain in agreement 
with experiments [8]. At the present stage, the 0(a 2 ) corrections to the other non-leptonic 
mode, b — > ccs, becomes one of the main issues with the theoretical predictions of the semi- 
leptonic branching ratio. The 0(a s ) correction to the c — s current in the mode b — > ccs is a 
monotonically increasing function of the mass of the charm quark, and the ratio of the 0(a s ) 
correction to the tree-level rate diverges in the threshold region 5 = (M&— 2m c )/2Mt, -C 1 due 
to the logarithmic term (In 5). Furthermore, although the fully analytic form of the 0(a 2 ) 
corrections, including the charm mass dependence, has not been obtained, it is known that 
the 0(a 2 ) corrections also contain the leading logarithmic terms (In 5) 2 and the next-to- 
leading logarithms (In S). We conjecture that the large corrections come from the threshold 
region and think it meaningful to elucidate the physical origins of the large corrections. 
With this goal in mind, in the present paper, we focus on the 0(a s ) and 0(a 2 ) corrections 
to the c — s current in the threshold region of the decay b — > ccs and estimate the quantum 
corrections which come from the two physical origins, the Coulomb force between the c and 
s quarks, and the anomalous dimension of the c — s current. 

It is well known that the Coulomb force between the proton and the electron produces the 
additional corrections to the tree-level decay rate of the neutron beta decay n — > pe~v e (see, 
e.g, the textbook [9]). The wave function of the electron is distorted by the Coulomb force 
generated by the electromagnetic charge of the approximately static proton. The distorted 
wave function forms the Fermi function in the decay rate. Considering the Fermi function 
in the neutron beta decay, here we attempt to incorporate the effects of the Coulomb-like 
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gluon exchanges between the c and s quarks in the threshold region of the decay b — > ccs. 
In Ref . [in] , it is pointed out that the Fermi function for the c — s current produces a next- 
to-leading logarithm in the 0(a 2 ) corrections. However, in that work, the Fermi function 
is not introduced to the decay rate formula which contains the phase space integrals of the 
final states. We believe it is worthwhile to obtain an estimate of the effect of the Coulomb 
distortion of the strange quark in the decay rate formula in order to obtain the more precise 
predictions of the Coulomb corrections. Considering this point, the main purpose of the 
present paper is to incorporate the Fermi function into the decay b —>■ ccs and to obtain 
the corrections which are produced by the Coulomb force between the c and s quarks in the 
threshold region. We then compare the Coulomb corrections with the Abelian parts of the 
perturbative corrections at the 0(a s ) and 0(a 2 ) levels. The other origin of the logarithmic 
terms is known. It is shown that the Wilson coefficients for the anomalous dimension of the 
heavy-light current, which is called the hybrid anomalous dimension, reproduce the leading 
logarithmic terms in the perturbative corrections of the two-body decay into one heavy and 
one light particle [TTJ [121 02] • Furthermore, when the Wilson coefficients are improved by 
the non-logarithmic terms at 0(a s ) in the perturbative corrections, the improved coefficient 
can produce the next-to-leading logarithmic terms in the 0(a 2 ) corrections [M f [T5l \16\ [17]. 
In the present paper, we confirm that the leading and next-to-leading logarithms in the 
perturbative corrections to the decay rate T(b — > ccs) can also be reduced to the Wilson 
coefficients for the hybrid anomalous dimension. 

The remainder of the present paper is organized as follows. In £j2]we obtain the threshold 
expansion of the perturbative corrections at 0(a s ) and 0(a 2 ) in the decay b — > ccs. In £J3] 
we obtain the Coulomb force corrections and compare them with the Abelian parts of the 
perturbative corrections. In £j4] we show that the Wilson coefficient for the hybrid anomalous 
dimension reproduces the leading and next-to-leading logarithmic terms in the perturbative 
corrections. In §|5] we give a summary. 



2 Results of the perturbative calculation 



We start with the effective Lagrangian for the decay b — > ccs 




C = -V cb V*-^$ c r(l - TS)^lWl - 7b)^c]. (1) 



2 The effective Lagrangian should include the Wilson coefficient which comes from the evolution of the 
renormalization group from the scale of the W boson mass to that of the bottom quark mass IH [5j [7] . 
But in order to compare the perturbative corrections with the Coulomb corrections, which is done in <|3l we 
do not need the Wilson coefficient, and we therefore omit it for simplicity hereafter. 
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We first calculate the tree-level decay rate by using the ordinal decay rate formula 

1 f d?p c 1 d 3 pc 1 d 3 p s 1 



T(b — ► ccs) 



2M h J (2tt) 3 2p° (2vr) 3 2p° 5 (2vr) 3 2p> 
(27r) 4 ^ 4 )(p 6 -p c -p g 



Pa)- 



(2) 



where .M is the invariant matrix and Mj, is the mass of the bottom quark. The momenta 
of the bottom, charm, anti-charm, and strange quarks are written Pb,p c ,Pc, and p s . The 
tree-level decay rate T(b — > ccs)q is obtained as 



ccs Jo 



Vl -4e 2 (l - 14e 2 - 2e 4 - 12e 6 



+24e 4 (l - e 4 )ln 



1 + VI - 4e 2 



.i-vT34?;j' (3) 

where the ratio of the charm quark mass to the bottom quark mass is written e = m c /Mb. 
The decay rate at e = is denoted by T M = N c \V cb \ 2 \V cs \ 2 G 2 F Mj^ /(192tt 3 ), with N c = 3. We 
can also obtain the decay rate from the so-called 'factorization formula', 



T(6 



CCS 



\V c b\ 2 \ V CS \ 2 G 2 F 



(M b -m c ) 2 



dq 2 
(27T) 



Here, X^ p and are defined as 



r(w~ -> cs) 



(4) 



(5) 
(6) 



where q and e^(q) are the momentum and polarization vector of the W gauge boson. When 
we calculate X^ v and Y^ v in Eqs. flSJ) and (jBj), we use the Lagrangians £ = t/> c 7 m (1— 75)^6 
and £ = ■0 s 7 At (l — ■y5)tp c W~ , respectively, where the coupling constants are normalized to 
unity. The mass squared of the boson, q 2 , is the variable of integration in Eq 
it runs over the region in which both of the decay processes b — > cW~ and W~ 
occur. We can calculate X^ v at tree level as 

dil 
An 



Ej), and 

cs can 



X 



[IV 



1 y/\(Mg,m*,q* 



47T 



M 3 



g^M 2 + g^{p b ■ q) - ie^q a p wp 



(7) 



where the variable Q represents the direction of the momentum of the W boson, and A is 
defined as A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + be + ca). We can also calculate Y^ v at tree level 
as 



AnJ~q< 



[A {q 2 g^ -q^)-B^q 



(8) 
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where An and Bq are the transversal and longitudinal parts, respectively, in the W boson 
decay. They are calculated as 



An — 



-fl-- 

£_f e 2 - 2 



1 + 



2u< 



where u 
form : 



= yff/Mb. Substituting Eqs. © and flU into Eq. 
r(6 — > ccs) 



(9) 
(10) 

we obtain the following 



6r, 



dw 2 v / Apy) 

u 2 {l + e 2 -u 2 )A + ({l 



and (1101) into Eq. fllip . we can obtain a result which coincides 
Thus, we can reproduce the tree-level decay rate in Eq. ([3]) from 
JJ. Since we focus on the threshold region of the decay b — > ccs 
in the present paper, we obtain the threshold expansion of the tree-level decay rate in Eq. 
(ED as 



Then, substituting Eqs.([£ 
with that given in Eq. ([3]) 
the factorization formula 



_,. . 4096 r7/2 / 261 . 1283 

T(b -> ccs) = T u 5 7/2 1 5 H 5 

v 1 u ^ 35 V 66 264 



231 
66~ 



1283 



(12) 



35 V 66 264 
where the expansion parameter 5 is defined as 5 = 1/2 — e. 

Next, we consider the estimation of the QCD corrections to the decay b — > ccs. In the 
present paper, we concentrate on the QCD corrections to the part of the decay W~ — > cs in 
the factorization formula (j4j)|f| We write the decay rate including the quantum corrections 
to the c — s current as 

2 



r(&^«*)=r (i + — Ax+f^) A; 

V 7T V 7T / 



(13) 



where Ai and A2 are the 0(a s ) and 0(a 2 ) corrections. We write the 0(a s ) and 0(a 2 ) 
corrections to the transversal part A Q in Eq. ([8]) as Ai and A 2 and denote the sum of the 
transversal part as A = A + A\ + A 2 . In the same way we write the longitudinal part as 
B = Bq + B\ + B 2 . The 0(a s ) corrections A\ and B\ are given as the functions of z = u 2 / e 2 
as [6] 

2z 



A 1 



C F A 



AC*) + 



2^+1 



h{z) 



B i = (^CfBoIMz)-!}, 



3 It is known that penguin contributions to the process b — > ccs exist and that they are small 
we ignore them for simplicity. 



(14) 
(15) 

Here 
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Figure 1: The corrections Ai and A2. In the plot of Ai, the solid curve and the dashed curve 
represent the exact results and the first three terms in Eq. (1221) in the threshold expansion, 
respectively. In the plot of A 2 , the solid curve represents the first three terms in Eq. 



where Cp = 4/3 and fi(z) and f2(z) are given by 

fi(z) = H + 2Li (1) + \nz ln-^ - ? ln(z - 1) + 



In h - ln(z - 1) 

z — 1 z 



z-l 
5 1 

2~ ~z 



Inz, 



1 /z-1 
+ 



z - 1 



biz. 



(16) 
(17) 



Here, Li(x) = — J ln(l — t)dt/t is the standard dilogarithm function. Replacing Ao and 
So with A\ and JBi, respectively, in Eq. (flT|) . we evaluated the integrals numerically and 
obtained the form of Ai plotted as the solid curve in Fig. [TJ From this graph, we can 
read off the value Ai = 4.46 at the reference point e = 0.3, which produces the corrections 
(a s /7r)Ai = 0.28 when we use the value a s = 0.2. This enhancement by about 30% is also 
reported in Ref. |6J . In order to observe the 0(a s ) corrections in the threshold region, we 
use the expansion parameter p = 1 — 1/z and expand A± and B\ in Eqs. (114"]) and (|T5|) up 
to the first three terms as 



A 1 



Bi 



71 



3L f 



-7T 



V 4 2 f 



+ P 



1 , 

-7T' 

3 



139 11 r 2 2 

1 L 7T 2 

18 3 p 9 



■3p 



23 
18 



3 p 
35 



2 ' 



(19) 



These expansions coincide with the results of the perturbative calculations presented in 
Ref. [19] . Here, we define L p = lnp. We change the variable of integration from u 2 to 



u 



[1/2 — u) in Eq. ( fill to obtain the threshold expansion. With A\ and B\ in Eqs. ([I 
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and (1T91) . we obtain the following threshold expansion : 



6r, 



J du(l - 2u) Vu + 6y/(l -6 + u)(l + 6-u)(2-6-u) 



2~ U 



[1 -5 + U + 5 2 -u 2 )A 1 + 



- + 5-5 2 
4 



u 

2 



r, 



7t 



387008 16384 



--5 + 5 2 
4 

n 16384 



2 

8192 



(20) 



1155 



315 



-7T 



35 



ln2-^L,U 7/2 + 



125825888 8192 



-it 



45045 45 

/ 22449032 2627584 „ 
f + 7T 2 



15015 



10395 



35 

1040384, n 520192 T \ rq/2 

^T5- ln2 + ^15- L T + 
727040^ _ 363520 x 5U/2 



99 



99 



+ 



(21) 



Here, we define L$ = ln<5. Dividing the quantity in Eq. (1211 by that in Eq. (fl2"|) . we obtain 
the threshold expansion for Ai as 

/6047 4 



+ -it 
2112 9 



4 ln2 -2L S ) + 



41184 

11236399 640 , 320 r \ r2 

1 ln2H L,5U 2 + - 

226512 99 99 / 



570307 128 , „ 
H ln2 



(22) 



We plot the first three terms in Eq. ( 1221) as the dashed curve in Fig. [TJ We can see in the plot 
that although the first three terms in Eq. (1221) are sufficient to reproduce the approximate 
value in the region e > 0.45, they are not sufficient at the realistic reference point e = 0.3. 

Finally, we obtain the threshold expansion of the 0(a 2 ) corrections. Their analytic forms, 
A 2 and B 2 , are not known. The first three terms in the threshold expansions of A 2 and B 2 
are obtained as [19] 



A 2 
B 2 



(v) C F -p 2 (C F A v A + C A A v NA + T R N L Al + T R A v H ), 
(^) 2 CV(1 - P)P 2 (C F A S A + C A A S NA + T R N L A S L + T R A S H ), 



(23) 
(24) 



where the explicit forms of A^, A^, and so on are given in Ref.[19] . Here, we define the 
coefficients as Cp = 4/3, C A = 3 and T R = 1/2, and we count the number of light quarks 
as N L = 3 in the present b quark decay. Using Eqs. (ITT]) . (Tl2|) . (1231) and §M§, we obtain the 
threshold expansion for A 2 as 



A, 



3367549 



+ 



42877 2 2 4 

7T 7T 

12672 81 



444217. 



-In2 



— vr 2 ln2 + 26(ln2) 2 

76032 12672 81 6336 27 v ; 

444217, 91 2 , , T 13 r2 70 > , ,i 

■Li - — n 2 L s + 261n2 L 5 + — L 2 - —((3) + 



12672 



27 
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1471797053 41175 2 7578347, 464 2 , 1088 „ x2 

-vr 2 + ln2 + vr 2 ln2 (ln2) 2 + 



8895744 9152 28512 243 9 

5 + 



7578347 T 160 2r 1088, „ r 272 r2 104 >/oN 
" L 5 + tttttVt 2 L 5 - — — ln2L 5 - —Lj + — C(3) 



57024 243 9 9 27 

29171145503 202553411 „ 10217309, 5552 2 , 

-vr 2 + ln2 + vr 2 ln2 + 



3669494400 12231648 141570 891 

60608,, „ x2 10217309 T 2656 2r 60608, „ r 

(In2) 2 H L 5 H it L$ H ln2 L a + 

891 v ; 283140 891 891 

15152 2 200 i 2 

-89r^ + -33-C(3)j«5 +•••• (25) 

We plot the first three terms in Eq. (125 j) in Fig. [TJ We can see in the plot that the divergences 
in the threshold region are more rapid than the divergences of Ai in Eq. fl22l . due to the 
existence of the terms L 2 . Although the corrections in Eq. (|25|) at the reference point e = 0.3 
are not valid, we read it as A 2 = 25.2 for order estimation. Using this value, the decay rate 
is estimated as T(b —>■ ccs) = T (l + 0.28 + 0.10), where the second and third terms are 
the 0(a s ) and 0(a 2 ) corrections in Eq. (fl3|) with a s = 0.2. Then, in order to make the 
magnitudes of the corrections in the threshold region concrete, we take the values Ai = 13.3 
and A 2 = 324 at the reference point e = 0.49. We estimate the corrections at the reference 
point as T(b ccs) = T (l + 0.85 + 1.31). This estimation shows that the sum of the first 
few terms in the perturbation series at the threshold cannot produce a valid approximation. 



3 Coulomb force correction 

In this section we obtain the corrections which are produced by the Coulomb force between 
the anti-charm and strange quarks in the final state of the decay b — > ccs. Our method to 
obtain the Coulomb force corrections is to replace the plane wave function of the strange 
quark by the wave function in the presence of the Coulomb potential. This method is 
used to incorporate the effect of the Coulomb distortion of the electron plane wave by the 
electromagnetic charge of the proton in the neutron beta decay rate [9]. The decay rate 
formula in Eq. ([2]) is not suitable for the distorted wave function by the Coulomb potential 
because the distorted wave function is not a momentum eigenfunction and we cannot extract 
the momentum conservation law which is expressed as the delta functions in Eq. 00 We 
choose a normalization of one particle per volume V for the plane waves in this section. 

4 If we can construct a wave packet which coincides with the Coulomb-distorted wave function near the 
origin and a plane wave far from the origin, we may be able to include the Coulomb force corrections in a 
more precise way. But in this paper, we do not try to construct such wave packet for simplicity. 
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Then we start the inclusion of the Coulomb corrections with the decay rate formula as 

n n ,g 

dT = 27r\(f\H int (t = 0)|i>| 2 S(M - n 7^ V - ( 26 ) 



i=l 



Here, the interaction Hamiltonian is defined as H^(t) = J d 3 xH^(x,t), where is 
interaction Hamiltonian density and M is the mass of the initial particle. Since the initial 
and final states of the decays are energy eigenstates, we can factor out the delta function 
which gives the energy conservation law. Using the interaction Lagrangian in Eq. (Tj[|) for 
the present decay, we can write the transition matrix as 

(ccs\H- mt (t = 0)\b) = V cb V; s ^= J d 3 x[^(l - 75 )M^(i - 75)#]- 

(27) 

In the ordinal perturbation theory, we use the plane waves for the participating particles as 



where the spinors of the particle and anti-particle are defined as 

u s {p) = ( £ s ) and v s (p) = (s) ( %f S ) , (29) 

with s — 1 for the up spin and s = — 1 for the down spin. Here we employ the standard 
Dirac representation, in which 7 is a diagonal matrix and the plane waves are normalized 
to one particle per volume V as J v d 3 x il) s (x)^ ip s ' (x) = 8 SjS i. In order to obtain the Coulomb 
corrections, we have the following two approximations. The first approximation is to ignore 
the small Pauli components {a ■ p)x s /{E + m) in the wave functions of the nonrelativistic 
charm and anti-cham quarks, because their small components do not contribute to the 
leading term of the threshold expansion. The second one is to replace the wave function of 
the strange quark in Eq. (|28|) with its value at the origin as 

r s '(x)^r s \x = 0), (30) 

because the contributions of the wave function at a radius far from the origin can be ignored 
in the first term of the threshold expansion. These two approximations are justified by the 
fact that the resultant decay rate coincides with the first term in the exact decay rate at 
tree level in Eq. f[T2"j) . as shown below. Taking advantage of the approximation in Eq. f[3"0"j) . 
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we replace the plane wave at the origin with the wave function in the Coulomb potential 
at the origin. We briefly review the solution of Dirac equation with the Coulomb potential 
(see, e.g, the textbook Ref. [20] ). Here we keep the finite mass of the strange quark unless 
stated otherwise. We solve the Dirac equation, 

(il^ - \e\QA^ - rnjtfj(x) = 0, (31) 

with the Coulomb potential Aq = — \e\Q s /r and A, = 0. To solve this equation we can use 
the general form of the wave function in a spherically symmetric potential, 



= ( f i!^i-) (32) 



where 3^' 3 j represents the spinor-spherical harmonics, and the total angular momentum 
is defined as j = 1+1/2, with the orbital angular momentum I. Although we have one 
more choice for the wave function, in which the orbital angular momentum of the upper 
two components is larger than that of the lower two components by unity, we use the wave 
function in Eq. (132|) . because the s-wave of the particle (strange quark) in the present case 
contributes to the Coulomb corrections. We obtain the solutions of the radial parts with the 
continuum energy spectrum as 



fi(r) 
9i{r) 



Nie- ikr (2kry- 1 x 

■r-\ -I I P. (fi -I- 1 1 -U On/ OohrW \ 

(33) 



V g+rg iPifc 1 + 2 7 , 2ikr) + j^iF^ + 1, 1 + 2 7 , 2ikr)} 
VW^fikFxte, 1 + 2 7 , 2ikr) - j^iFxit + 1, 1 + 2 7 , 2ikr)} 

where \F\ is the Kummer function and we have 7 = yJT+ l) 2 — cto, £ = j + iv, v = 
a^E/h, a = otQin/k and cto = — \e\ 2 QQ s . Using the normalization condition 

'/•'•'• ^ 3 (xM' J Hx) = h^'^Kk - k% (34) 
we fix the normalization constant as 

N _ e ,u,2 . |r(i + 7 + ^)| 

|r(i + 2 7 )| ' {6b) 

where T(z) is the Gamma function. We can rewrite the solutions of the radial parts in Eq. 
(!33j) near the origin as 

\9i{r = R) )~ Nie [2kR) \ VE=fH[l - j^] J ' (36) 
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where we define the radius near the origin as R ~ 0. When the Coulomb potential is switched 
off and the Coulomb wave is reduced to a free spherical wave, the inner product of the wave 
functions in Eq. (132]) with j = 1/2 near the origin is given by 

which differs from that of the plane waves, ^ s (xyi/) s (x) = 1/V. In order to replace the plane 
wave in Eq. ( 130]) with the Coulomb-distorted wave in Eq. ( 132]) near the origin in a consistent 
way, we modify the normalization constant of the Coulomb wave in Eq. (l35j) as 



For this purpose, here we define the Fermi function as 



F(E,a ) 



|^ 3 (x,«o)| 2 



\^ 3 (x,a 



(39) 

r=R 



with R ~ 0. Using Eqs. (|32l) . ( J35l) . and (|36l) . we obtain the explicit form of the Fermi 
function as 

F{E, a,) = 4^±^ . {2kR) ^-» . <~ . (40) 
E + m |r(l + 27o)p 

with the relations E = \J m 2 + k 2 and 70 = \/l — a^. Hereafter, we omit the mass of the 
strange quark again and write the energy of the strange quark as E s = \p s \. Substituting 
the plane waves in Eq. (128]) into Eq. fl27j) . except for the strange quark, and using the 
renormalized Coulomb wave in Eq. ( l38i) for the strange quark, we obtain the transition 
matrix element as 

(f\H int (t = 0)\i) = V cb V: s ^^{2*f5\p c + mry3\cs) 

Hx r '^ l X r )3\cs)} (41) 
= {f\H- mi \i)v + {f\H iTA \i) A , (42) 



where (f\H^\i)v and (f\H^\i)A are defined as the first and second terms in Eq. (jH 
referring to the vector and the axial vector parts of the bottom-charm current, respectively. 
We write the anti-charm and strange quark current as j^{cs) = 1 I r s 7 M L'0 g = $^7^0^, where 
the left-handed fields are written as g = Lip 5 and $ s = L\l/^-' 3new , with the projection 
operator L — (1 — 7s)/2. We average the square of Eq. ()42l) over the spins and obtain the 
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contributions of the vector and axial vector parts as 

2 1^ \\j\ H mt\ % )v\ ~ \Vd>Vca\ y2 y b {^s,a ), (43) 



r.r'.s,. 



2 lul-^intPMl = |KbV^| -pj y 3F{E s ,a ), (44) 



where the interference term of the vector and axial vector parts does not contribute to 
the spin-averaged decay rate. When we derive Eqs. (143j) and (1441) . we use the formulae 
£>g(aO^(aO + - L / 2 ^ and £ s $ s (f)t$ s (f)| r=R = F(£7 a ,a )/K Substituting Eqs. gHD 
and (1441) into Eq. (|26|) . we obtain the decay rate including the Fermi function F(b — > ccs)° 
as 

r s 

r(6 -> ccs) c = 768 / d^ 2 (l - 2z)V5 - z - 5 2 + z 2 ■ F(E S , a ), (45) 

Jo 

where the variable of integration is defined as z = E s /2Mb. We compute only the contribu- 
tion to first order in S and obtain the decay rate as 

/1019 7T 2 



r(6 -> ccs) c = F c 



1 + vra + {-^q - le + y - HM b R5) - 41n2 )a, 



(46) 



Here, we have used the expansion of the Fermi function in the coupling constant, F(E, «o) = 
1 + 7rao + (tt 2 /3 + 3 — 7 e — \\i{2ER))al + • ■ •, and written the tree-level rate as Tq = 
r M (4096/35) 8 7 ^ 2 , which coincides with the first term in Eq. (1T21) . Here, 7 e ~ 0.577 is Euler's 
constant. Since the Coulomb corrections originate in the effects of the Coulomb-like gluon 
exchanges between the c and s quarks, we can expect that the 0{olq) and 0(olq) corrections 
in Eq. (|46p are included in the Abelian-type one- and two-loop diagrams shown in Fig. [2J 
In the correspondence, we can introduce the strong coupling constant as a = Cpa s and 
rewrite the series in Eq. (T46l) as 

F(b — > ccs) c _ _ o(ot s \ / 1019 7r 2 , „ „\ 

~r§~ = i + ^(-) + ( w - 7e + T -MA W -4h l2 ) 

x Cl*<(^)\ (47) 

Here we compare the corrections in Eq. (T4T1) with the abelian part of the perturbative 
corrections in Eqs. (122]) and (1251) . In the comparison we assume that the correspondence 
of terms is categorized according to the powers of the quantities, n, C(3), and logarithmic 
function. For example, the 7r 2 Ls term at 0(a 2 ) in Eq. (1471) is compared with the n 2 Ls term 
in the abelian part of the perturbative corrections in Eq. ( 1251) . rather than with the L$ term. 
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Figure 2: The Abelian-type diagrams which include the Coulomb-like gluon exchanges 
between the anti-charm and strange quarks. 



The coefficient Cfk 2 in the 0(a s ) correction is three times larger than the similar tt 2 term, 
47r 2 /9, in Eq. ( 122]) . We conjecture that this difference is compensated for other corrections 
which have the other physical origins. At (9(a 2 ), we extract the Abelian part from A2 in 
Eq. (1251) by taking only A^ and A S A in Eqs. (1231) and (124"1) and thereby obtain the first term 
as 



The term —Cp7r 2 Ls in the 0(a 2 ) corrections in Eq. (J57J) accounts for 120% of the similar 
7i 2 Ls term, — (5/6)Cpii 2 Ls, in Eq. (148]) . Because the L 2 and L$ terms are singular at the 
threshold, these terms are not cancelled by other non- logarithmic terms. Based on the 
assumption that the terms are categorized according to L$, 7r 2 , and so on, we conjecture 
that the sum of the Coulomb corrections and other tt 2 Ls corrections, which come from other 
physical origins in the Abelian theory, is — (5/6)Cpir 2 Ls in Eq. (1481) . The numerical impact 
of the Coulomb corrections, however, is not clear at the moment. The differences among the 
other terms, 10197r 2 /210, 7r 4 /3, and — 47r 2 ln2, at the 0(a 2 ) in Eq. (14TI) and the corresponding 
terms in Eq. (148]) are large. Although the Coulomb corrections do not contain terms of the 
forms ln2, (ln2) 2 , L$, ln2L,5, L 2 , and C(3), the perturbative corrections do contain such terms. 
Further, the Coulomb corrections contain j e , while the perturbative corrections do not. We 
also believe that these differences are compensated for the other corrections. 

4 Hybrid anomalous dimension 

In this section, we show that the leading and next-to-leading logarithmic terms in the first 
order of the threshold expansions in Eqs. (1221) and ( 1251) originate in the anomalous dimension 



abel 




(48) 
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of the heavy-light c — s current j = qTQ in the heavy quark effective field theory (HQET), 
which is called the hybrid anomalous dimension [IT] . Here, q is a massless quark field, Q is 
a heavy quark field in HQET, which satisfies the relation Q = 7°Qj and T represents the 
Dirac gamma matrices. First, solving the renormalization group equation 



(^-7/»(a.(p)))c(^,M) = 0, 



(49) 



we obtain the Wilson coefficient for the hybrid anomalous dimension as 



C(n, M) = C(M, M) exp 



("Or 



(50) 



where the hybrid anomalous dimensions in HQET are given in Refs. [Til [I~2l [13j [TT] and [21 



as 



with 7^ 



(o) 



-4 and % 



(i) 



^= + ^-ni — Mir 2 . The beta function are given by 



(51) 







4vJ 



(52) 



with fto = 11 — (2/3)n/ and /5i = 102 — (38/3)n;. The hybrid anomalous dimension in 
HQET does not depend on the structure of the gamma matrices T. We write the initial 
condition in Eq. ( 1501) as C(M, M), which depends on the structure T. In order to deal with 
the universal part in Eq. (1501) . we set the initial condition as C(M,M) = 1 hereafter. Valid 
initial conditions are introduced later through the matching coefficients. Using Eqs. fl5Tl) 
and ([521), we obtain C(/i, M) in Eq. (EDI) as 



C(E s ,m c 





2/3 

exp 


r tT, 


'if 


/3A 


-a s (m c ). 




V7f 


ft/ 



(53) 



where we write fi as the threshold momentum E s , like the energy of the strange quark in 
the presently considered decay, b — > ccs, and M as a heavy quark mass m c , like the mass of 
the charm quark. Furthermore, we substitute a s (E s ) to Eq. fl53|) with the running coupling 
constant up to two-loop order, 



a s (m c ) 



1 + a s (m c )fa log ^ + a s (m c ) 2 #^ W ^ ' 



(54) 
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and obtain the square of the Wilson coefficient in Eq. (153]) as the following series in powers 
of the strong coupling constant a s (m c ) : 



C 2 (E s ,m c ) 



1 + 



a s (m c ) 



TT 



An * 



L 2 



7T 



+ 



(55) 



The leading logarithmic terms in the C?(a s ) and C^cti?) corrections in Eq. (I55l) with = 3 
coincide with those in the perturbative corrections in Eq. (|22l) and (l25l) when 5 is replaced 
with E s /m c . This agreement shows that the leading logarithmic terms in the perturbative 
corrections come from the hybrid anomalous dimension. By contrast, the next-to- leading 
logarithmic term at 0(a 2 ) does not coincide with that in Eq. (j25|) . This disagreement 
between the next-to-leading logarithmic terms is natural for the following two reasons. The 
first reason is that the next-to-leading logarithms in the transversal part, A/A Q , and the 
longitudinal part, B/Bq, in the factorization formula (llip are different, while the leading 
logarithmic terms are the same. The second reason is that although the product of the 
leading logarithmic term and the non- logarithmic terms at 0(a s ) can produce the next-to- 
leading logarithmic terms at 0(a 2 ), the square of the Wilson coefficient in Eq. (!55l) does 
not contain non-logarithmic terms at 0(a s ). In connection with the first of these reasons, 
the relations between A/A , B / B and the universal part, R(p), are given by [17] 



MP) 
B{p) 
Bo(p) 



C v (m c yR(p) 
mJrric 



C s (m c ) R(p) 



where the matching coefficients are obtained as [22] 



C v (m c ) 

C s (m c ) 
rhJiric) 



TT 



C F f a s {m c ) 
2 V TT 
1 _ ^( a s( m cj 
3 



(56) 
(57) 

(58) 
(59) 
(60) 



m c 6 v TT 

Using Eqs. (Q and (ITS]) and the leading and next-to-leading logarithmic terms in (123]) . we 
can extract the universal part R(p) in Eq. (I56j) as 

17 4 Q i /«,(m> 2 



R(p) 



1 + 



a s im c 



TT 



-2L + 



-TT 



p ' 3 ' 9 
1657 97 



+ 



36 



TT 

— TT 2 + Til (— + —TT 2 ) ) L 

27 V 6 27 1 1 p 



(61) 
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which is consistent with Eq. (l57j) . In connection to the second reason stated above, regarding 
the lack of non-logarithmic terms, we consider the non-logarithmic terms at 0(a s ) in Eq. 
fl6"Tl) and introduce the following factor at the scale of the threshold momentum E s : 



The relation between a s (E s ) and a s (m c ) appearing here is given in Eq. (l54"j) . When the 
coefficient C 2 (E s ,m c ) is improved by the factor D(E S ), it coincides with the universal part 
in Eq. (16T1) . up to the next-to-leading logarithms at 0(a 2 ), as§ 



This means that the next-to-leading logarithmic terms in Eq. (1251) can be reduced to the 
hybrid anomalous dimension and the running coupling constant in Eq. (I54"l) . 

5 Summary 

We have studied the physical origins of the corrections to the c — s current of the decay 
b — > ccs at the threshold. We first obtained the 0(a s ) and 0(a 2 ) perturbative corrections 
in Eqs. (I22p and (I25p . respectively. Second, we obtained the corrections in Eq. (14 7p . which 
are produced by the Coulomb force between the c and s quarks. The Coulomb corrections, 
Cpir 2 at 0(a s ) and — Cp7c 2 \n 5 at 0{a 2 ) 1 account for 300% and 120% of the terms with the 
same forms in the Abelian parts of the perturbative corrections, respectively. The differences 
among the other terms in the 0(a 2 ) corrections in Eq. ( )47l) and the perturbative corrections 
of the same form are large. These differences should be compensated by other corrections 
which have other physical origins. We finally confirmed that the Wilson coefficient for the 
hybrid anomalous dimension reproduces the leading logarithmic terms in the perturbative 
corrections. We also confirmed that when the Wilson coefficient is improved by the non- 
logarithmic terms in the 0(a s ) perturbative corrections, the improved coefficient can produce 
the next-to-leading logarithms in the 0(a 2 ) perturbative corrections. We have identified 
three challenging problems in the present paper. The first problem is to reproduce the 
Coulomb corrections in Eq. f|4Tl) with perturbative calculations. In such calculation we 
should calculate only the loop diagrams shown in Fig [2] and extract only the soft-gluon 
contributions in the loop integrals. The second problem is to determine the accuracy of the 

5 In Ref.[17] , The universal part, R(p), is obtained by using Eq. ([49]) and the 0(a s ) corrections of R(p). 
In the same way, here we obtained the relation (I63|) by introducing the factor in Eq. (|62|) which has the 
non-logarithmic terms at 0{a s ) in R{p) and is defined at the energy scale E s . 




(62) 




(63) 
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method used to incorporate the Coulomb force corrections in £j3j Solving this problem, we 
will be able to identify those terms that come from the Coulomb force corrections with the 
higher precision. The last problem is to find the relations between the two next-to-leading 
logarithmic terms of the Coulomb corrections in Eq. (|4"T|) and the Wilson coefficient for 
the hybrid anomalous dimension in Eq. (ISTl) . We should observe only the Abelian part of 
the hybrid anomalous dimension. Although the two next-to-leading logarithmic terms may 
originate in the soft and hard regions of the loop momentum, respectively, these relations 
are not trivial. We need more sophisticated studies to reveal these relations. 
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